Equivalence between discrete quantum walk models in arbitrary topologies 
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Coin and scattering are the two major formulations for discrete quantum walks models, each 
believed to have its own advantages in different applications. Although they are related in some 
cases, it was an open question their equivalence in arbitrary topologies. Here we present a general 
construction for the two models for any graph and also for position dependent transition amplitudes. 
We then prove constructively their unitary equivalence. Defining appropriate projector operators, 
we moreover show how to obtain the probabilities for one model from the evolution of the other. 

PACS numbers: 03.67.Lx, 05-40.Fb 



I. INTRODUCTION 

In essence, there are three different possible implemen- 
tations for quantum walks (QW), all taking place in dis- 
crete spaces (graphs). Two are quantum analogues of 
Markov chains [H,T2] : generalizations of diffusion- like dy- 
namics, where time is continuous (CTQW) [3[, and dis- 
crete time unitary maps [l|, known as coin QW (CQW). 
Recently, it has been shown to exist a direct relation be- 
tween these two cases [3, H|. The third (SQW), also a 
discrete time formulation, is physically appealing since it 
is based on the idea of scattering in multi-port interfer- 
ometers 

Due to the importance of Markov chains in random 
algorithms, CQW and CTQW are extensively used in 
their quantum versions j8j, often being more efficient 
because the exponentially faster hitting times of QW 
[3, 0, E3] ■ Relevant is also the finding that one can imple- 
ment universal quantum computation through scattering 
processes in a CTQW model . Regarding applications 
for SQW, they seem to be particularly suitable [12J to 
solve searching problems [13j ]. 

It is generally believed that CQW, CTQW and SOW 
have distinct advantages in different contexts @, EH, Ell • 
For instance, universal port gates in quantum computa- 
tion like Hadamard's [161 ] are easier to implemented with 
CQW, the most common formulation in the mathemati- 
cal and computer science literature. On the other hand, 
the construction of SQW in graphs of arbitrary topolo- 
gies is more direct, and thus conceivably simpler to real- 
ize experimentally. Furthermore, since it is based on an 
scattering approach, analytical techniques are far more 
developed for SQW than for coin models fl4j . 

Hence; (i) to formulate in a constructive way both 
the CQW and SQW for graphs of any topolo gy [17j and 
for position dependent quantum amplitudes [la |: (ii) to 
prove their unitary equivalence in the general case (pre- 
viously established only for particular situations @, [T^]); 
and finally (iii) to show how to obtain the probabili- 
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ties for one model from the other; are fundamental be- 
cause the following. First, a constructive rather than an 
abstract formulation makes easier concrete implementa- 
tions of such systems, even for arbitrary topologies (e.g., 
by means of BE condensates in optical lattices [20J ) . Sec- 
ond, such results would bridge the gap between two ap- 
parent distinct quantization schemes for a same class of 
systems. Finally, they also would show that the usages 
for one model are equally possible for the other. The 
points (i)-(iii) are then the goals of the present paper. 



II. GRAPH STRUCTURES AND BASIC 
DEFINITIONS 

We assume an undirected simple arbitrary graph [13], 
whose nodes are labeled in Z. Its topology is entirely 
determined by the sets Vj = {ji,32, ■ ■ ■ )3nAi which rep- 
resent the Nj nodes connected to the node j. Thus, if 
ji belongs to Vj, then there exists exactly one edge be- 
tween j and ji. Also, to any node j we associate the set 
of integers Aj — {1,2, . . . , Nj}. Each element a of Aj 
corresponds to a different edge attached to j, in a one- 
to-one relation. Note that if j r and j s have a common 
edge, then there are two integers numbers, a r and a s , 
associated to such edge, one due to j r and other to j s . 

A given mapping (function) on a graph is said locally- 
adaptable if: (i) it can be constructed for any graph node 
j; (ii) for each j, it depends only on the edge structure of 
the nodes in Vj ; and (iii) it is always well-defined regard- 
less the number of elements in Ay j . This is an important 
concept because if one can establish time evolution re- 
laying only on locally-adaptable mappings, then the re- 
sulting dynamics is valid for any graph topology. 

So, consider two locally-adaptable functions that direct 
reflect the specific structure of a given graph. The first, 
e : Aj — > Vj, associates each a from Aj to a single ji 
from Vj, such that e(a;j) gives the node connected to j 
through the edge labeled a with respect to j. For the 
particular example schematically depicted in Fig. 1 (a), 
we have Vj = {31,32,33,34} and Aj = {1,2,3,4}, thus 
e(>i;j) = ji, e(cr 2 ;j) = 32, e(cr 3 ; j) = h, and e(cr 4 ;j) = 
ji, where each Oj assumes one of the values in Aj. The 
second, 7 : Aj —> Ay., maps each (Ji from Aj to the 
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FIG. 1: Examples of locally-adaptable mappings, (a) Map- 
pings which are directly associated to the graph topology, (b) 
Arbitrary mappings, used to define the CQW evolution. 



single element a of A e ( CT( .j) representing the edge joining 
j and e{<Ji]j). Again in Fig. 1 (a) we have 7(01; j) = at, 
7O2; j) = cr;, 7(03;.?') = cr m , and -y(a 4 ; j) = a n . Here 
crfc belongs to A JX = {1,2}, cr; to A J2 = {1,2,3}, and so 
forth. Note that in general 7(7(17; j); e(cr;j)) = cr. 

The above mappings are strictly related to the 
graph specific topology. However, more general locally- 
adaptable functions can also be defined. The next three 
will be very useful. /1 : Aj — * A\> j extends 7 since it 
associates each a from Aj to an unique arbitrary ele- 




FIG. 2: In the coin (a) and scattering (b) QW formulations, 
a is associated, respectively, to nodes and edge states, (c) An 
example of relabeling for scattering states, Eq. (JTJl. 

ment of A e ( CT; j), e.g., in the particular case of Fig. 1 
(b), we have fx(ar,j) — <Jh, where ah is in Aj 1 = {1,2}. 
Although such function is locally-adaptable, for an ap- 
propriate and consistent latter construction of the quan- 
tum evolution along the whole graph, we consider an 
extra restriction for /1. From its very definition, for 
any a G Aj we have that fi(j(a;j);e(a;j)) is also 
an element of Aj. Then, we impose additionally that 
M(7(ov; j);e(oy; j)) ^ ^(a s ;j);e(a s ;j)) if a r ^ a s 
(a r , a s in Aj), i.e., the set {^(a; j); e(a; j))} = Aj. 
Observe that this restriction can always be fullfiled what- 
ever the graph topology. For instance, in Fig. 1 
(b) we have {/x(cr n ; j 4 ), ^{a k ; ji), /i(er/; j 2 ), n(a m ; j 3 )} = 
{<?i> <?2> 03 j C4} = Aj. Naturally \l induces two other 
locally-adaptable functions. Indeed, suppose ay run- 
ning over Aj, so Vtj = {{a t ,ji)} = {{^{a r ;j),e{a r ;j))} 
(at e Aj i and ji e Vj, (a = 1, 2, ... , Nj)) is a set where 
to each pair (ai,ji) corresponds a distinct a G Aj. We 
have thus v : f2j — * Aj and a : Clj — > j, such that 
v { (J i\ji) — a an d a<{&i]ji) — j- By construction 

j = a(n{a; j); e(a; j)) = e(u(a;j);a(a;j)), 

a = v(fi(a;j);e(a;j)) = /j,(v(a;j);a(a;j)). (1) 



III. THE TWO DISCRETE TIME 
FORMULATIONS 

For the coin version, the states are defined on the graph 
nodes j [2l[. Thus, the cr's labeling the edges attached 
to j can be associated to the quantum numbers repre- 
senting the different "outgoing" directions leaving j, as 
schematically shown in Fig. 2 (a). Hence, we have as the 
base states {|j, cr) c }, where for each j, a = 1,2, ... , Nj 
and (a' ', a") c — 8j> j>> cv a /i , which spans the Hilbert 
space 7i — L 2 (Z x Zjv z ). 
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To establish the system dynamics, we first consider the 
shift operator S (and its adjoint [13] ), such that 



S\j,o) c 



\e((r;j),n(a;j)} c 



(2) 



From Eq. (JTJ) , it follows that 5* 5 = S = I in H. 
Then, for each j, let (7W to be a "coin" operator, repre- 
sented by a iVj x A^ unitary matrix, whose action over 
a basis state (of quantum number j) is C^\j,a) c = 



J2a> 3 =i c a> > a\j' a ') c - Finally, we set the unitary one step 
time evolution as [231 



,0') 



^ = s££<? (j) b>><<M| c , 



(3) 



know at time n what is the probability pj (n) to be in the 
position state j (which means a node (edge) in the coin 
(scattering) model), regardless the value of the coin (di- 
rection) quantum number a. So, we define the scattering 
and coin projector operators as 

P s U ' a) = \j,^){^j\ s + \e((r;j),j(a;j)}(j(a;j),e(a;j)\ s , 



JV, 



= £|j»<<7,j|c 



(5) 



<j=i 



The desired probability is thus the expected value 

p^(n) = (*(n)|P|*(n)), |*(n)> = C7"|*(0)) (6) 
for P one of the expressions in Eq. ([5]). 



which is valid for any topology (encoded in the functions 
e and 7) and defines a very general time evolution for the 
problem through the functions fj,, v and a. 

For the scattering version, note first that even for a 
same graph, the a labeling [24| for the scattering [f| 0] 
can be completely distinct than that for the coin formu- 
lation (see Figs. 2 (a) and (b)). Thus, in principle the 
functions e and 7 can assume different values in the two 
cases (and they will be distinguished when necessary). 

Now, two quantum states are defined along each edge, 
e.g., for the edge connecting the nodes j and j m in Fig. 
2 (b), we denote the state "incoming" to the node j 
(jm) by |j, a s ) s (|j m ,0i) s ), which can be written also 
as |e(cr i ;j m ),7(cr i ;j m )) s (|e(<r s ; j), 7(cr s ; j)) s ). Actually, 
for any edge, if one state is \j,a) s , the other is given by 
\e(a;j)^(a;j)) s . So, the basis set is {|j,<t)} s , spanning 
the Hilbert space TL — L 2 (Z x Zjy z ), as in the coin case. 

For the dynamics, we set U s = R + T [23[, where the 
action of the operators R and T are given by [22j 



£ 4 J ?Je(a;i),7(«;i)) s , 

a£ Aj , a^cr 

r ^ji^;j) l e <>; i)> 7(0-; j))s, 



|e(<r;i),«) s .(4) 



We also define AT, x A, scattering matrices rw, such 

(for both a' ^ a in 



that T, 



and ri J ?„ = t U) 



a' ,a 



Aj). If for all j we impose that rw is unitary, then the 
coefficients r and i satisfy the usual relations in scattering 
theory [25l. |26|. So, U s is unitary. 



IV. OBTAINING THE PROBABILITIES 

For QW, stochasticity - in the classical sense - comes 
into play only through measurements, when one calcu- 
lates the probabilities for the walker to be found in dif- 
ferent locations along the graph [27| , Suppose we shall 



PROVING THE EQUIVALENCE OF THE 
TWO FORMULATIONS 



For so, three steps are necessary: (a) to establish a 
correspondence between the different walks states; (b) 
to properly associate their time evolutions; and (c) to 
construct projector operators to obtain the probabilities 
of one in terms of the other. 

Regarding (a), note that we always can define a locally- 
adaptable function ip : Aj — > Aj, which for each node j, 
maps the quantum number a associated to a specific edge 
in the scattering formulation to the quantum number o' 
labeling the same edge, but in the coin formulation. For 
example, for the situation in Fig. 2, we have ip(ak'ij) — 
cti, <p(cri;j) = ct 2 , <p(cr r ;j) = cr 3 , and ip(a s ;j) = ct 4 . Also, 
the actual cr's values are not relevant. They are just 
a way to label nodes and edges states. Thus, without 
loss of generality, for any j we always can rename one of 
the model states by |j, er) => \j, 4>(&] j)) for <t> a bijection 
Aj — > Aj. Choosing to retag the scattering case, we 



consider the following particular <f> 

<P{v,j) = v(lc(<p(<r;j);j);e.c((p{(T;j);j)), 



(7) 



whose "action" is pictorially represented in Fig. 2 (c). 

Then, using this new notation for the scattering states, 
one has the isomorphic unitary operator E : TL — > TL ((J 



E \j,°)s = \j,o) c - 



(8) 



associating the scattering state a incoming to j to the 
coin state a outgoing from j (see Fig. 2 (c)). 

For point (b), from the explicit form of U c and U s , one 
finds that by setting one of the two akin relations 



p« 



C U) = r 0') 



(9) 



then, in both models the time evolution transition prob- 
ability amplitudes are exactly the same. In this case, the 
resulting dynamics are unitary equivalent since 



U s = £t u E . 



(10) 
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Finally, for (c) even if the two formulations are unitar- 
ily related, the resulting probabilities - through projec- 
tions - are not 6]. This is so because each description 
assumes distinct spatial configurations, nodes or edges, to 
characterize the system, thus not leading to same prob- 
abilities. In fact, in each edge the two scattering states 
are mapped by E to coin states in different nodes. Thus, 
the probability, Eq. ©, to be in a unique node is not 
equal to the probability to be in a unique edge. 

However, there is a very direct way to obtain the walk 
probabilities for the coin (scattering) model from the 
scattering (coin) model. We just define 

P U)\ =Etp<J)E, P r W)| =EPWE*. (11) 

S|c C ' C|g S \ / 

Then, suppose we construct a SQW with arbitrary r's 
and i's. By using P, '^ of Eq. © into Eq. ©, we get 

the scattering walker probabilities at the step n. But now 
(i) i 

if we use P s | for this system, the resulting probabilities 
are exactly those from a coin model, for which the coin 
matrices elements are given by such r's and i's values 
according to the correspondence in Eqs. ©. The other 
way around, to get the scattering model results from the 
CQW, follows in the same fashion. 

VI. REMARKS AND CONCLUSION 

Although there are some few discussions in the lit- 
erature on how to formulate discrete random walks in 
general terms [HI, E3, EH , here we have developed a ex- 



plicit procedure that allows one to write, in a construc- 
tive way, the time evolution operator directly from the 
system topology and local dynamics. It also includes the 
case of position dependent quantum amplitudes (through 
the C^'s and T^'s). Moreover, in the CQW case, it is 
not necessary all the matrices C^> to have the dimensions 
equal to the largest coordination number of the graph, as 
in certain formulations [141 ] . 

In the present framework, a regular graph can be de- 
fined in terms of the A^'s sizes (e.g., all equal to N) and 
the features of e, 7 and fi (e.g., independent on the j's and 
having specific patterns along the graph). For instance, 
for all nodes with the same number of edges, we know 
that for the coin, Ti. can be written as the direct product 
of two subspaces, thus \j,a) — > \a) ® \ j). So, we natu- 
rally find from Eq. © that £i E^l C^\j,a)(<r,j\ = 
EjCWQUWlaa in 0. 

Usually, it is believed that analytical methods are eas- 
ier to implement for the SQW than for CQW 0, HI]. 
Since in fact they can be mapped each other, the exist- 
ing methods for the former should be extendable to the 
latter. Our constructive approach may serve as a guide 
to implement such extensions [30| . 

Lastly, by using the cross operators in Eq. (fTTj) . we 
have been able to calculate - for different examples - the 
probabilities for SQW and CQW from a single implemen- 
tation. This will be communicated elsewhere [301 ] . 
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